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Abstract 

We conduct extensive numerical experiments aimed at finding the 
admissible range of the ratios of the first three eigenvalues of a planar 
Dirichlet Laplacian. The results improve the previously known theo- 
retical estimates of M Ashbaugh and R Benguria. We also prove some 
properties of a maximizer of the ratio A3/A1. 



1 Introduction 

Let be a bounded domain in W 1 , n > 2. We consider the eigenvalue 
problem for the Dirichlet Laplacian, 

-Ait = Xu in 0, (1.1) 



u\an = 0. (1-2) 

Let us denote the eigenvalues by Ai(f2), A2(f2), • • • , (we will sometimes 
omit explicit dependence on f2 when speaking about a generic domain), 
where < Ai < A2 < A3 < . . . . The corresponding orthonormal basis of 
real eigenfunctions will be denoted {uj}'jL 1 . 

For the last fifty years, the problem of obtaining a priori estimates of the 
eigenvalues and their ratios has attracted a substantial attention. The ex- 
isting results can be roughly divided into two groups — universal estimates, 
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valid, as the name suggests, for all eigenvalues and all the domains in R n , 
which do not take into account any geometric information, and isoperimetric 
estimates for low eigenvalues. We briefly survey some known results below; 



the reader is referred to the very detailed survey paper |Ashf and references 
therein for a full discussion. 



Universal estimates 

Probably the first, and best known, estimate of this type is the Payne-Polya- 



Weinberger inequality [PPW]. 



A m +1 < A m + 



mn 



(PPW) 



This was subsequently improved by Hile and Protter [HiPr], and, in 1990's, 
by Hong Cang Yang [Yan], whose implicit estimate 



m+l 



\m+l 



H — I A- 

n 



< 0. 



(HCY) 



remains the best universal estimate so far for the eigenvalues of the Dirichlet 
Laplacian. 

The general method of obtaining ( [PPW ) and (HCY), as well as of similar 
estimates for a variety of other operators, has been the use of variational 



principles with some ingenious choices of trial functions, see [Ash]. Recently, 
an alternative abstract scheme, based on the so-called commutator trace 
identities, which easily implies, in particular, ( PPW| ) and ( HCY ), has been 
developed in [LevPai], see also [HarStu]. 



By their very nature, the universal estimates are generically non-sharp. 



Isoperimetric estimates 

Both ( pPWp and ( pCY[ ) give, for m = 1, the estimate 



Ai n 

This upper bound cannot, in fact, be attained. Already, Payne, Polya, and 
Weinberger conjectured that the actual optimal upper bound on the ratio 
of the first two eigenvalues of the Dirichlet Laplacian is 



Ai Ai 



i 2 

■Ai/2,1 



A 4 

n-dimensional ball Jn/2— 1,1 



for Q C 



(AB 



(here j p>q stands for the q-th. zero of the Bessel function J p (p), so, in a planar 



case n = 2, K 2 « 2.5387 compared with (|PPW| ) bound -^(0) 

Ai 



ncR 2 



<3). 
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Conjecture ( ABp ) was eventually proved, only in the early 1990s, by 
Ashbaugh and Benguria [|AshBenl| , |AshBen3|1 , using, in particular, sym- 
metrization techniques going back to the Faber-Krahn inequality, 



Ai(0) > Ai(n*), 



where fi* is an n-dimensional ball of the same volume as f2. 

We would like to mention, in this context, extensive computational ex- 
periments designed to verify ( |ABq ) by Haeberly [ Hae , HaeOveH . 



Statement of the problem 



As mentioned, ( [ABqD gives the full description of the range of the possible 

A 2 

values of ratio of the first two eigenvalues of the Dirichlet Laplacian, — 

Ai 

A 2 

for domains in the Euclidean space (the obvious lower bound is — > 1). 

Ai 

In fact, similar results were also obtained for domains in S n and M n . A 
natural extension would be to find optimal upper bounds on the range of 

(A2 A3 
— , — 
Ai Ai ( 

in particular for planar domains. In other words, we would like to find, for 

x := — - and v := -A, the function 
Ai Xi 



y*{x) :-- 



max 



£(n) 



(1.3) 



and the number 



Y* ■- 



max y*(x) 
xe[l,K 2 ] 



max -r~(P) 
Ai 



(1.4) 



or their best possible estimates. We will use notation ( |1.3| ) and ( |1.4| ) when 
looking for maxima in particular classes of domains as well. 

Despite an apparent simplicity of this problem, and a wide attention it 
has attracted, it turned out to be rather difficult. In | AshBen5 , [AshBcng ], 
Ashbaugh and Benguria proved a complicated upper bound for y*(x) and 
also demonstrated that 



35 

3.1818 ps — < Y* < 3.83103 . 



(1.5) 



Their estimates improve upon previous results due to Payne, Polya, and 
Weinberger, Brands, de Vries, Hile and Protter, Marcellini, Chiti, Hong 



Cang Yang, and themselves; see [ AshBen5 , AshBenC] and their earlier papers 



[AshBen2, AshBen4| for extensive bibliography and details of proofs. We 
present their estimates and other known facts in the next Section; just note 
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at the moment that the lower bound in is attained when $7 is the 

rectangle R a := [0, 1] x [0, a] with a = J^. 

In the current paper, we describe extensive numerical experiments aimed 
at improving ( |1.5|) , We also show, using perturbation techniques, that the 

rectangle R rj does not maximize the ratio — and indicate a class of do- 

V 3 ^1 

mains among which a possible maximizer could be found. 
Acknowledgements 

The original suggestion to conduct numerical experiments on the range of 

— — , — - J for planar domains came from Brian Davies; we would like to 
Ai Ai J 

thank him, as well as Mark Ashbaugh and Leonid Parnovski, for valuable 
discussions and advice. It is the authors, nevertheless, who take full re- 
sponsibility for the realization of this idea, and any criticism for possible 
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2 Known results for the range of ( — ^, — ^ J for pla- 



nar domains 
Explicit solutions 



Ai'Ai 



The spectral problem (LI), (L2) admits a full solution by separation of 
variables when Q, is, for example, a disjoint union of a number of rectangles 
or circles. For a reference, we collect below the results on the range of 

— -, — - J in these cases. 
Ai Ai / 

Rectangles Let R a := [0, 1] x [0, a] be a rectangle with the side ratio a; 
without loss of generality a > 1. Then 

a7 ™ - l^FTi 



and 



a n - for a > \l — 
a? + 1 ~ V 3 

4 ° 2 + 1 f 1 <r <r 
lor 1 < a < 



U 2 + l 

Thus, for rectangles, in notation (|1.3|) and ([L 



y(x) — y (x) | rec tangles 



( 8 5 20 
—x for 1 < x < — , 

11 - - 2 



4 



and the maximum value of ^~ is 

Ai 



Y*\ V) 



I rectangles 
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attained when a = \ —■ Note that for this particular rectangle A3 is a 
V 3 

degenerate eigenvalue: X^(R ij) = X^(R /—), and it is the only a for 

V 3 V 3 

which Xs(R a ) is not simple. 

In the (x, y)-plane, (2.1) corresponds to the two straight lines inter- 



/20 35 

secting at the point I — , — 

Circles For a single circle, x = y = K 2 . As easily checked, for a union of 
more than one disjoint circles of arbitrary radii, 

y*(*)lcircles^2, for 1 < X < K 2 . (2.2) 

Its graph in the (x, y)-plane is a straight line parallel to the x-axis. 

Disjoint unions The following easily checked fact shows that one cannot 
obtain higher values of y*{x) by considering disjoint unions of sets 
from two different classes. Namely, let Cj be two arbitrary classes of 
domains, with corresponding functions y*(x)\ c . (not necessarily de- 
fined for all x £ [1,-^2])- Then, for any domain £1 = S7i U ^2 with 

^2 ^3 

Qj G Cj, we have, for x = — (O) and y = —($7), the inequality 

Ai Ai 

y < max (y*(x)\ Ci , y*(x)\ C2 , K 2 ) . 

Other domains There are other domains, like sectors of the annuli, el- 
lipses, etc., for which the problem of funding the eigenvalues is reduced 
by separation of variables to the problem of solving some transcenden- 
tal equations. However, the latter one is often not easier than the 
numerical solution of the original problem, so we do not treat these 
cases here. 

The graphs of y(x) = y*(x)| rectangles and y*(x)| circles are shown in Fig- 
ure ||. 

Ashbaugh-Benguria estimates 



In [AshBenf ], Ashbaugh and Benguria proved, using a wide variety of meth- 



ods, the following upper bounds for y*(x): 

y*{x)<K 2 x for 1< x < 1.396", (ABi 
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/*(x) < 1 + x + ^2x - (1 + x 2 )/2 for 1.396" < x < 1.634", (AB 2 ) 



y*(x)<F(x) for 1.634" < x < 1.676", 



< H(x) - x for 1.676" < x < 2.198" 



and 



y*(x)<G(x) for 2.198+ < a; < 2.539" 
where the functions H(x),F(x) and G(x) are defined by 
6 for x = 1 



(AB 3 ) 
(AB 4 ) 

(AB 5 ) 



. / 4/3(/? + 7 ) 2 (x - l)(x - fhftf + 7 ~ I)) 2 

i<TA, V V (2/3 - 1)(2 7 - l)(x - - 0(4x - 2 - r? - £) 
for x > 1 



with /3 = 7/ + \/rj 2 — r] and 7 = ^ + y^ 2 — £, ^(x) is the middle root of the 
cubic 

2xy 3 - 2(5x 2 + 3x + l)y 2 + (6x 3 + 39x 2 + 2x - l)y - (24x 3 + llx 2 - 4x - 1) = 
and 



G{x) = inf 



ft 1 



+ 



1-/37(2/3-1) 



/3>i/ 2 V2/3-l C 2 (/3)(x-/?7(2/3-l))-l, 
with the infimum taken over values of /3 satisfying a; > (5 2 /{2f3 — l) + l/C2(/3) 



and with 



J0.1 



C 2 (/3) 



2/3-1 



t 3 J 2(3 (t)dt 



/3 io.i 

/ Ul p {t)dt 



We remind that Jo(^) denotes the standard Bessel function of order zero 
and j'0.1 is its first positive zero. For the derivation of the bound ( IAB4I ) and 
more discussion of it, see [AshBen5]. The other bounds given above are due 
to Hong-Cang Yang (see | Yan | for ( AB2D ) and Ashbaugh and Benguria (see 
[|AshBen2| 1, [|AshBen4| for ( |ABl| ), |AshBen5| for ( [AB^l) , and [|AshBen6|l for 
(IAB3D). 



The admissible region for 



A2 A3 
Ai'Ai 



defined by ( ABq| ), obvious bounds 



2 > 1 and — ^ > and inequalities ( ABj )-( AB5 ), is shown in Figure ||. 



Ai Ai Ai 

We make two remarks following | AshBenq| : 
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4.2 



Ashbauch-Benguria estimate 

I Admissible region 

4 - — Rectangles 

Circles 




for comparison are the maximum values of — as functions of — for rectangles and 

1 1 

disjoint unions of circles. 



Remark 2.1. The inequalities ( ABi[ )-(A.Bs) all apply on broader intervals of 
x-values than the intervals specified explicitly with them; the given intervals 
indicate the range for which the corresponding inequality gives the best 
bound yet found. 



Remark 2.2. The absolute maximum of the right-hand sides of ( ABi )— (SB 



occurs at the point where F(x) has a maximum within the interval where 
it is the best bound. That happens at the point (x,y) « (1.65728,3.83103), 

and imply the best upper bound (1.5) yet proven for 



A3 
Ai ' 



3 Numerical analysis of random domains 

To the best of our knowledge, there have been no large scale numerical 
experiments on low eigenvalues of the Dirichlet Laplacian for planar domain. 

(A2 A3 
— j — 
Ai Ai 

we conducted such experiments for a variety of domain classes. 
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General method 



For each particular domain, the calculation of the first three eigenvalues has 
been conducted using a standard finite element method implementation via 
PDEToolbox iPDEToolj and FEMLAB jFEMLABfl in Matlab, with two or 

three mesh refinements. For simple domains with relatively "high" values 

A 3 

of the ratio — , optimization with respect to the parameters describing do- 
Ai 

mains of the particular class was performed in order to maximize this ratio. 
The results of calculations for some classes of domains are described below 
and summarized at the end of this Section. 

For each class of the domains we represent the results in the following 

A 2 . 

graphical form. The range [1,-^2] of possible values of x = — is split into 

Ai 

subintervals of length 5x (normally ~ 0.05). In each subinterval we choose, 

if it exists, a domain with maximal y = — ^ and plot the corresponding 

Ai 

point (x,y). For comparison, the graphs of y{x) = y*(^)| rec tangles and/or 

V*( X ) (circles are sllOWn - 



Triangles, quadrilaterals and ellipses 

We start with cyclic calculations for all triangles, with angle step 2.5°. For 
the triangles with relatively high ratio of A3/A1 we repeat the procedure 
in the local neighbourhood with angle step 0.5°. The results are shown in 
Fig. § 

The computational procedure for quadrilaterals is essentially the same 
as the one for triangles, with parameters a, (3, 7 and 5 in the region (0, 7r) 
(see Fig. ||). We choose an angle step of 2.5°. Note that quadrilaterals with 
negative a and /3, or 7 and 5 do not have to be considered separately - 
they fit into the scheme above by choosing another diagonal as a starting 
point and re-scaling. In case of relatively high ratio A3/A1 (> 3) we repeated 
the calculation with an angle step 0.5° in the local neighbourhood of that 
quadrilateral. The results are shown in Fig. |3[ 

The results of cyclic calculations for the ellipses, with axis ratio varying 
between 1 and 5 with step 0.1 (0.001 in the vicinity of the ellipse with highest 
A3/A1), are also shown in Fig. ||. 

Remark 3.1. Rather surprisingly, Fig. || suggests that 

V ( X ) (quadrilaterals ~ V ( X ) (rectangles (3-1) 

We give a partial explanation of this fact in the next Section, see Remark |4.6| . 
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Figure 2: Parametrization of a quadrilateral. 
Annuli and random sectors of annuli 

The calculations for annuli with with inner radius 1 and outer radius r 

A3 

demonstrated that the corresponding value — (r) is monotonically increas- 

Ai 

ing from 1 to K2 as r changes from 1 to 00 (although convergence, for large 
r, is very slow — just logarithmic). These results are not very informative 
and we do not include them in the graphs or the summary table below. 

In calculations for sectors of the annuli of angle 9, we choose r randomly 
in the interval (1,20) and 8 randomly in the interval (0.0l7r, 1.997r). The 
results of calculations are shown in Fig. || 

Pseudo-random polygons 

For polygons with more than four vertices, cyclic calculations through all 
possible values of the geometric parameters with some reasonable step be- 
come impractical due to the time constraints. Instead, we choose to perform 
calculations for randomly generated polygons. We employ the following 
simple procedure for generating a pseudo-random polygon with ./V vertices 
vi, . . . , vat lying inside a square [0, l] 2 . 

Vertices Vi,V2,V3 are chosen randomly using any pseudo-random genera- 
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Figure 3: y*(x) for triangles, quadrilaterals and ellipses 



tor. 

Vertices v,-, j = 4, . . . , N — 1. We choose a possible vertex at random. If 
the interval [vj i; Vj] intersects any of previously constructed sides 
[vfc_x, Vjfc], k = 1, . . . , j — 2, then we make another random choice. 

Vertex v^r is constructed in the same manner, but we additionally check 
that the interval [vjv, vi] does not intersect any of the existing sides. 

To avoid infinite loops, we abort the construction if the number of attempts 
at some stage exceeds some sufficiently big number (say, 200). We also put 
in place a restriction forbidding very small angles (which require special 
efforts in mesh generation). 

The collated results of calculations for pseudo-random pentagons, hex- 
agons and decagons are shown in Fig. ^. These results also include experi- 
ments on random perturbations of the rectangles constructed in the following 
way: N points were randomly chosen on the sides of the rectangle R a , with 
a G (1,5), and 1 < N < 8, and these points and the four vertices of the 
original rectangle were randomly moved by a distance not exceeding 0.1a to 
form an (N + 4)-gon. 
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Figure 4: y*(x) for random sectors of annuli and pseudo-random polygons 

Star-shaped domains (simply and non-simply connected) 

The procedure described above for the polygons does not work very effec- 
tively for polygons with large number of vertices — it often takes a long 
time to generate a suitable vertex Vj with j > 10. Thus, in these cases we 
restrict ourselves to star-shaped polygonal domains which are much easier 
to construct. Namely, for the vertices v = re 1 ® , we choose the angles ran- 
domly between and 2tt, and the radii r randomly between given numbers 
r± and T2- We conducted a series of experiments with a fixed number of 
vertices (13, 17 and 23), as well as a series of runs where the number of 
vertices was chosen randomly between four and thirty. 

Additionally, we conducted a series of experiments of non-simply con- 
nected domains of the types R \ S\, S2 \ R rj, and S\ \ S2, where Sj are 

V 3 V 3 

random star-shaped polygons such that Si C R rj C S2, and R rj is the 

V 3 V 3 

rectangle with the maximum — — . 

The results for pseudo-random star-shaped domains are collated in Fig. pi 
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Figure 5: y*(x) for pseudo-random star-shaped domains 



Dumbbells and Jigsaw pieces 

By a dumbbell we understand a domain of the type 

([0, 1] x [-h, h]) U C((0, 0), n) U C((Z, 0), r 2 ) , (3.2) 

where l,h,ri,rz are positive parameters, and C(v, r) denotes a circle with 
radius r centred at v. By a jigsaw piece we understand a domain of the 
type R\C, where R is a rectangle, and C is a circle with a centre "near" 
the boundary of the rectangle. Typical dumbbell and jigsaw piece domains 
are shown in Fig. |(| 

The results of numerical experiments on dumbbells and jigsaw pieces, 
with cyclical/random choice of the parameters, is shown in Fig. |7[ For 
dumbbells, we also optimized over the parameters for domains with A3/A1 ~ 
3.2, allowing additionally the centres of the circles to move in the vertical 
direction along the sides of rectangles. However, this did not lead to the 
improvement of the results. 

Summary of the numerical experiments 

We summarize the results of our numerical experiments in the following 
table. 
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Figure 6: Typical dumbbell and jigsaw piece domains 




Dumbbells 
Jigsaw pieces 
Rectangles 
Circles 



1 1.2 1.4 1.6 



2.6 



Figure 7: y*(x) for pseudo-random dumbbells and jigsaw pieces 



As seen in the last column of Table |], in each class of domains the 
A3 

maximum of the ratio — is attained, within the accuracy of computations, 

Ai 

on a domain with degenerate eigenvalue A3 ~ A4. The same, of course, holds 
for rectangles, see This allows us to conjecture that the absolute 

maximum and any local maxima of A3 ~ A4 are also attained on domains 
with degenerate A3. We give a partial proof of this conjecture in the next 
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Type of domains 


No. of experiments 


Y* 


h 


Triangles 


2145 


2.827 


0.016 


Quadrilaterals 


13222 


3.183 


1.6 • 10~ 5 


Sectors 


360 


3.149 


4.6 • lQ- b 


Ellipses 


142 


3.167 


9.1 • lQ- b 


Random polygons 


26867 


3.189 


6.1 • 10~ 4 


Random star-shaped domains 


18320 


3.159 


0.022 


Dumbbells 


2871 


3.202 


1.1 • 10" 4 


Jigsaw pieces 


1420 


3.178 


0.003 


Total 


65337 


3.202 


1.1 • 10~ 4 



Table 1: Summary statistics for numerical experiments. The value in the fourth 

column is 5^ := — ^— — -(O*), where f2* is the domain which maximizes the ratio 
A3 

A3 

— — = Y* in the corresponding class of domains 
Ai 

Section. 

The computed absolute maximum ratio Y* ~ 3.202 is attained on the 
dumbbell-shaped domain (|3.2j) with I = 1, h = 1.4510, r\ = 0.7814, and 
T2 = 0.7818, see Fig. ^. Note that the maximum value Y* is only slightly 
higher than the corresponding value ^*| rec tangies ~ 3.1818. 

Remark 3.2. Additional experiments were conducted in order to check whether 
a maximizer is likely to be a simply connected domain. Namely, for the 
dumbbell-shaped domain Q described above, we computed the eigenvalues 
for a number of domains obtained by removing a small hole from Q. In 
all the cases, the ratio of the third and the first eigenvalue for a perturbed 
problem was quite significantly less than that for Q. 

The graph of the function y*(x) built on the basis of all numerical ex- 
periments is shown in Fig. |||. 
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- - Rectangles 
Circles 



'1 1.2 1.4 1.6 „ 1.8 2 2.2 2.4 2.6 



Figure 8: y*(x) for all computed domains 



15 



4 Asymptotic results 



In this Section, using standard perturbation techniques, we establish several 
results which, although don't give the full answer to the question of maxi- 
A 3 

mizing the ratio — (0) among all planar domains f2, give some indication 
Ai 

which of the domains may or may not be a maximizer. We first proof the 
following 

A3 

Theorem 4.1. The rectangle R does not maximize the — among all 

V 3 Ai 

planar domains^ 



This should be compared, however, with Remark 4.6 below. 
We also give a proof of the following more general result, which justifies 
the remark made at the end of last Section. 

A3 

Theorem 4.2. Suppose that S!CR 2 is a local maximizer of — — (S7) among 

Ai 

planar domains with sufficiently smooth boundary. Then As($l) = \i{Q). 
We should emphasize here that neither the statement not the proof 



(found below) of Theorem 4.2 is fully rigorous. In the former, we do not dis- 
cuss the requirements on the smoothness of the boundary and the concept 
of a local maximizer; we also do not prove that maximizers actually exist. 
In the latter, we rely on the following unproven, although very plausible, 
conjecture. 

Conjecture 4.3. Let A3 be a simple eigenvalue of the Dirichlet Laplacian 
on a planar connected domain. Then not all nodal lines of the corresponding 
eigenfunction are closed. 

Such a conjecture is not unreasonable since, in general, it is quite difficult 
to construct domains for which even one nodal line of a low eigenfunction 
is closed, see [ |HoHoNa| 1 . 



Before giving the proofs of Theorems |4.l| and [4,2| , we recall, without 



proof, some classical results from the domain perturbation theory. The 



details can be found, e.g., in [ELel, SHSP]. 
Domain perturbations 

For simplicity, we restrict ourselves to domains in R 2 ; all the results stated 
here hold in any dimension. 



* We were informed by Niculae Mandrache that he had independently obtained a similar 
result. 
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Consider, for small values of real parameter |e|, a family of bounded 
domains Q £ in R 2 of variable x = (x\,X2), which are transformed by the 
change of coordinates 

x = x + eS(x) (4.1) 

into the domain f2 = f2° in R 2 of variable x. We assume that the boundary 
d£l and the vector-function S are sufficiently smooth. 
Let n be the outer unit normal to <9f2, and denote 

/ = S-n 

(in fact, ef is a smooth function on d£l which gives, up to the leading order 
for small e, the normal distance between dtt and d£l £ ). 

Denote by Ai < A2 < • • • < Xj < • • • the eigenvalues of the Dirichlet 
Laplacian on f2 and by {uj} the corresponding basis of normalized orthog- 
onal eigenfunctions (which are chosen real). Also, denote by A^ the eigen- 
values of the Dirichlet Laplacian on Q £ . For sufficiently small |e|, the A^ are 
continuous functions of e and tend to Xj as e — > 0. 

The following two results go back to Rellich. 

Proposition 4.4. Let Xj, j > 1, be a simple eigenvalue of the Dirichlet 
Laplacian on Q. Then X s - has the asymptotic expansion 

X £ j = Xj + eX jt i + e 2 X jt2 + . . . , (4.2) 

as e — ► ; where 

hi = -J f 

an 

The situation is slightly more complicated when Xj = ■ ■ ■ = Aj+ m is an 
eigenvalue of multiplicity m + 1. For simplicity, we consider just the case 
m = 1. 

Proposition 4.5. Let A& = A^+i be a double eigenvalue of the Dirichlet 
Laplacian on £1. Then, as e — > 0, A| and A| +1 still have the asymptotic 
expansions ( [4,2j ) (j = k, k + 1) with 

Afc,i = - min(e^i,e^ 2 ) , Afc+i,i = - max(e / ui, e^ 2 ) , (4.4) 
where jjL\, ^2 <ire two real roots of the quadratic equation 

(F k , k + v)(F k+1 , k+1 + fj.) — Fl k+1 = (4.5) 

and 

an 



duj 



dn 



da . 



(4.3) 
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A 7 

We will be in fact interested in the asymptotic expansion of — which 
_ A i 
follows from (O): 



A, 



+ 



Af Ax (Ai) 



rCAj-iAi-AriAjO + OCe 2 ), 



(4.7) 



Proof of Theorem 4.1 



Let Q = R /-g- be the rectangle ^ (a;i,a;2) : < x\ < 1, < X2 < \j )■ We 



shall construct an explicit perturbation S7 £ using (4.1) such that the first 
correction term in the asymptotic formula (4.7) is positive for e > 0, and 

therefore -|- 

Let 



> — 1 for sufficiently small positive e. 
Ai 



ST = I (xi, x 2 ) : < xx < 1, < z 2 < W - + eg(x\) 



where 



g{x\) = cq + v^c; cos(7r/xi) 



z=o 



We will choose the coefficients q later. 

The corresponding function / appearing in the asymptotic formulae 
above is 



f(xi,x 2 ) 




, < xi < 1 



if X2 



if (xi, x 2 ) G 8ft , X2^ \ %■ 



Note that we shall use (4.3) for computing A^i and (4.4) for computing 
X^i and A^i, since A3 = A4 is a double eigenvalue of the unperturbed prob- 
lem. Elementary but tedious calculations show that the correction terms 
Afc,i, k = 1,2, 3, 4, depend only upon the parameters Cj with j = 0, . . . , 4. 
For brevity, we omit the explicit expressions. 

Let us choose the parameters cq, . . . , C4 in such a way that Xs ; i = A^i 
(i.e., A| remains a double eigenvalue up to the linear terms in e). This, by 
Proposition 4.5, happens when ^3 3 = i<4 4 and ^3 4 = 0, which in turn leads 
to the following conditions on coefficients cy. 



C3 = ci , 



Under conditions (|4.8|) , asymptotic formula 
and becomes 

A| A3 _ 96^ 
Af ~ AT ~ 121 



c 4 = 9 C2 - 8V2c . (4.8) 
]) simplifies dramatically, 



(02 - \/2co)£ + 0(e 2 ) , 
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and we can choose cq and C2 in such a way that its right-hand side is positive 



for sufficiently small positive e. This proves Theorem iA 

Remark 4-6- Let = R a be any rectangle and consider the perturbations 
Q £ as above but with function / linear in x\, x% (and, naturally, a replacing 

8 throughout). Thus, we are considering quadrilaterals il e which are 
small perturbations of the rectangle R a . The same elementary calculations 



A 2 



then imply that, for x £ := — (J7 e 

Ai 



A 



— we obtain, up to and 
Ai 



and y e : 

inclusive of the terms of order e, that 

rectangles ' 

where y* (^) I rec tangies * s gi ven by the right-hand side of (|2.1|) . In other words, 
up to the terms of order e the rectangles are local maximizers among all 



quadrilaterals which are sufficently "close" to them, cf. Remark I3T 



Proof of Theorem 4.2 



Suppose that is a planar domain with sufficiently smooth boundary which 
locally maximizes the ratio ^ in the following sense: for any sufficiently 
smooth perturbation O, 6 determined by (4.1) we have 



A| < A3 



Ai 



(4.9) 



Assume additionally that A3 is a simple eigenvalue of the Dirichlet Lapla- 
cian in the unperturbed domain £1. We shall show that this assumption leads 
to the contradiction with Conjecture |4,3| , 

Since both Ai and A3 are simple eigenvalues, the asymptotic formula 



(4.7) becomes, in accordance with Proposition 4.4 , 



A, 



A3 
Ai 



(Ai 



/ A 3 



dn 



du?, 



dn 



da + 0(e 2 



Now, as e can be chosen both positive and negative, (4.9) can hold only if 



/ A 3 



m 



dn 



du 3 



dn 



da = 0, 



and since / is an arbitrary smooth function, this requires 



A3 



dn 



Ai 



dm 



dn 



everywhere on d£l. But the normal derivative of the first eigenfunction of 
the Dirichlet Laplacian is non-zero everywhere on the boundary, so the last 
formula implies that the third eigenfunction has the same property, and 



therefore all its nodal lines are closed, in contradiction with Conjecture 4.3. 
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5 Final Remarks 



On the basis of the numerical computations and the results proven above 
we make the following conjecture, most of which is still to be established (or 
disproved) rigorously: 

A3 

The domain maximizing the ratio — for planar domains is close in shape 

1 [71 

to the optimal computed dumbbell-shaped domain shown in Fig. is simply- 
connected, and has a smooth boundary. The maximal admissible value Y* 
is approximately equal to or is slightly greater than 3.202. 




Figure 9: Domain maximizing — on the basis of computations 
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